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We present a method for generating new deformed solutions starting from systems of two real 
scalar fields for which defect solutions and orbits are known. The procedure generalizes the approach 
introduced in a previous work [Phys. Rev. D 66, 101701 (R) (2002)], in which it is shown how 
to construct new models altogether with its defect solutions, in terms of the original model and 
^ ' solutions. As an illustration, we work out an explicit example in detail. 
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^ . PACS numbers: ll.10.Lm, 11.27. +d 

I. INTRODUCTION 

Kinks, domain wall, vortices, strings and monopoles are all well known examples of defect solutions with 
topological profiles for field theories in different dimensions [3, 0, Q ■ These solutions have been intensively 
studied since the seventies in high energy physics. The interest has been continuously renewed in different 
branches of physics since defect solutions in general appear in models of condensed matter as well as string 
theory. In particular, systems of real scalar fields have attracted attention with very distinct motivation, 
since they can be used to describe domain walls in supergravity [J] and braneworld scenarios with an extra 
dimension [5[. Therefore, it is important to find explicit analytic solutions for this kind of systems and is in 
this direction that the method conceived in ref. [(J has shown to be very useful. In that work it was shown 
that knowing a defect solution of a scalar field model for a single real field is enough for generating an infinity 
of new models with its solutions, all written in terms of the original model and solutions. 
\ For field theories involving two real scalar fields, the mathematical problem concerning the integrability of 

the equations of motion is much harder, as one deals with a system of two coupled second order nonlinear 
ordinary differential equations, and the configuration space shows a distribution of minima that allows for 
T^j- . a number of topological sectors. One way of simplifying the problem is to consider potentials belonging to 
the (wide) class corresponding to the bosonic sector of supersymmetric theories. This kind of systems can 
C*"** , be studied, via the introduction of a superpotential, in a first order formalism which allows (stable) BPS 
configurations [?|]. Even in this case finding explicit solutions can be a highly non-trivial task and therefore, 
any method for obtaining new solutions would be of great utility. 

For models with two interacting components, the solutions on each topological sector determine orbits in 
the configuration space, which can be expressed as a constraint equation 0(xi,X2) = 0. Based on this fact it 
was introduced in *§] a procedure called trial orbit method consisting in shooting an orbit and testing it on the 
equation corresponding to the model considered. Later, this method was adapted [9( for the searching of BPS 
states of systems of first order ODEs, leading to some advances. In recent years, other general methods for 
the investigation of complicated nonlinear problems arising^ in many-field systems which comprise multidefect 
solutions have been developed; see, for instance, [T3, EH E2, EH EH EH and references therein. Models of two 
scalar fields have also been used to describe complex phenomena such as the entrapment of topological defects; 
see, for instance [H, E3> EH • 

As will be shown below, the field deformation method introduced in [6| for one-field models also works for 
connecting ODE systems of two first-order equations in models with two scalar fields. The equations arising in 
this extended procedure are in general much more complicated than their counterpart for a single scalar field, 
given that not one but two deformation functions are now required, and it happens to be difficult to realize 
which pair of deformation functions would do the job in the right way, generating a well behaved deformed 
potential and solutions consistent with the equations of motion of the system. To overcome this difficulty, we 
take into account the fact that the actual solutions connecting two vacua of a topological sector live restricted 
to orbits in field space. Therefore, by deforming the first order equations for a two-field system while imposing 
the orbit constraint, we assure the consistency of the solutions of the deformed model at the level of the 
dynamical equations. 
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The paper is organized as follows. In the next Sec. II we briefly review the deformation method for one 
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component systems, and we extend the deformation procedure to two-field models. In Sees. Ill and IV we 
show how the incorporation of orbit constraints allows to obtain consistent deformed solutions for two-field 
models. Then we work out an explicit example in detail to illustrate the procedure in Sec. V. We end this 
work with some comments and conclusions in Sec. IV. 



II. THE DEFORMATION METHOD 



We start with a system of a single real scalar field x(x M ) : R 1 ' — ► M, o a bi-dimensional Minkowski 
spacetime (/x = 0, 1 with x° = xq = t, x l = —X\ = x), described by a Lagrangian with the usual form 

c = \d^x-v(x), (1) 

where the potential V(x) specifies the model. 

For static configurations (x — x( x ) '■ ^ — > the equation of motion reads 

fx = dvbd m 

dx 2 d X ' 

and the energy functional associated to the static solutions is given by 
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E[ X ] = J dx 



1 (dx 

dx 



V( X ) 



(3) 



Requiring for the energy of the solutions to be finite results in the boundary conditions dx/dx — > and 
V(x) as i ±oo. Thus, the physical solutions are constant at infinity, and their asymptotic values are 
minima of the potential. Performing a first integral of the equation of motion under these conditions, we get 
to the first order equation 



dx^ 2 
dx 



mx)- (4) 



In this work we will restrict our study to solutions of this first order equation which present topological 
(kink-like) character, in the sense that they connect two different minima of the potential. 

It is convenient to consider potentials of the form V(x) = h[W (x)] 2 1 where the prime means derivative with 
respect to the argument, and the functional W{\) is the superpotential. This allows us to write equation (j4]) 
as the gradient flow equations of W 

d f = ±W'(x) (5) 
dx 

Let us now describe the deformation procedure for a single real scalar field model, as introduced in @ . The 
prescription is the following. First, we define the deformed potential as 

tu^ YUM ™ 

m) = WW' (6) 

where / is the deformation function. 

This new potential determines the model for the deformed field (j> through the deformed Lagrangian 



C = -d^d^-U(4>)- (7) 



The first-order ODE determining the defect profile in the deformed model reads 



=2U M> (8) 
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with U{<f>) given by ©. This can also be seen as the gradient flow equation of the deformed superpotential 
WO) defined by 



Second, we connect solutions of the original and deformed models through the deformation function 

X (x) = f[4>(x)\, (10) 
Here the deformation function f{4>) is assumed to be bijective - see however [Tlj |. 

Therefore, if defect solutions x{ x ) °f © are known, the link between the two models provides the defect 
solutions 4>{x) of ([8|) by just inverting the field transformation (fTO]) 

^) = /- 1 (xW)- (ii) 



III. EXTENSION TO TWO-FIELD MODELS 



At the level of the first order equations, the deformation method can be directly generalized to models with 
two scalar fields. We start with the model 

£=^x-d»x~V(x), (12) 

where x is an isospinorial real field x — (xii X2)- Suppose that the potential energy density can be written in 
the form 

, 1 fdW\ 2 1 fdw\ 2 
F( Xl , X2 ) = -(-] j (13) 

where the superpotential W is a well behaved function in the space of scalar fields x(x,t) £ Maps(R 1:1 , R 2 ). 
A subtle point is the following: Because (|13[) is a PDE equation there can be several independent solutions 
for W - not merely a global change of sign as in the one-field case - see Then the static finite energy 

solutions (topological defects) of this model satisfy the first-order ODE system 

(14) 





dW 


dx 


dxi 


dX2 


dW_ 


dx 


dX2 



Now we choose a deformation function / : R 2 — > R 2 such that 



Xi(x) = fi((j>t(x),<f>2(x)) 
X(x) = f(ct>(x)) { , (15) 

X2{x) = f 2 {4>l(x),(t>2(x)) 



where <j>i and (j>2 are the deformed fields. Then the first-order ODE system (fT4|) becomes 

dih d -t + d 2 h d -t = d h W 



(16) 



where 
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This system can be rewritten as 

d -t = lk (92f2d h W-d 2 hd h W) 



dx J(/) 

(17) 

d -t = 777) (d l f 1 d h W~d l hd h W) 

where J(f) — <9i/i d 2 f 2 — 82/1 d\f 2 (the Jacobian of /) , and we assume that J(f) 7^ 0. (Note that this fact 

can be relaxed by restricting the deformation method to act on a family of open sets in R 2 where J(f) 7^ 0). 
Equations (|17p can be interpreted as the first-order ODE system 

(18) 

d ~t = 9 2 W 

where W(4>i,(j)2) is the superpotential of the deformed system and the right hand side derivatives denote 
In terms of the original superpotential, the deformed one is determined by the PDE system 



{ jfc(d2f2d h W-d2fld fa W) = dxW 

V W) ( d ^^ d h w - d ihd fl W) = d 2 W 



(19) 



Therefore, if any two-component defect solution x(x) °f 03} i s known, the link between the two models 
provides two-component defect solutions <fi(x) of (fT5| by calculating 

Vice versa, if two-component defect solutions 4>(x) of (|18p are known, the link between both models provides 
two-component defect solutions x(x) of (|14[) by applying the transformation 

X(x) = f($(x)) 

The existence of f~ l is associated to the above mentioned precisions about the zeros of </(/) via the Inverse 
Function Theorem. 

At this point, the differences between working with one or more scalar fields appear. First, the ODE system 
(fT8t is supposed to give the solutions for the deformed system. This will be true only when those solutions 
satisfy, besides the first order system, the dynamical equations for the deformed potential 

V((h,<h) = ^PiW) 2 + (9 2 W) 2 ]. (20) 

This is automatically satisfied by any solution of (fT8|) . But beside this, one is assuming that the deformation 
leads to well-behaved, smooth potentials. Therefore, we should ask for the first order derivatives of the 
deformed superpotential W(4>i, </>2) to be continue or, alternatively, for its second order cross-derivatives to be 
identical c^W = d 2 iW. Imposing this condition on the systems (fT5|) and (fTT)]) leads to a very complicated 
constraint which suggest no obvious choice of the deformation functions fi and f 2 - For this reason, in order 
to make progress we need to introduce some assumptions to simplify the situation to a tractable case. We will 
consider functions of the form /1 = fi(4>i) and f 2 = f 2 (4> 2 ). This significantly reduces the complexity of the 
constraint, which becomes the simple condition 

{{d x hf - (d 2 f 2 ) 2 ) d hh Wdxh d 2 f 2 = (21) 

This is true whenever /1 and f 2 satisfy difi — ±d 2 f 2 but, as fi{4>\) and f 2 {4> 2 ) are functions of different 
fields, this expression seems to be nonsense. However, since the solutions 4>\{x) and 4>2{x) live on an orbit in 
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configuration space, condition (|2ip must be understood as a function of the solutions, therefore we rewrite it 
as 



dfr(fa) _ ± df 2 (fa) 



(22) 

orbit 



When the deformation functions depend on one field only, the PDE system (flU]) takes the simpler form 

(23) 

2 W(fa,fa) = 7T-? ' 

C2/2 

Each one of these equations resembles the deformation recipe applied to the case of a single scalar field model, 
so one could think that it is just a duplication of the one-field procedure. However, the superpotential W now 
depends on both deformation functions fi(fa-) and f^fa), as a consequence of the interaction between the 
fields, and then the deformed superpotential W depends on both fa and fa, and the resulting deformed model 
describes interacting fields. 

As we have shown, requiring that the deformed model has a smooth potential leads to a condition on the 
deformation functions. Using this, and the fact that the fields are enforced to obey an orbit we are able to 
construct a consistent deformation procedure, detailed in the following section. 

IV. THE ORBIT-BASED DEFORMATION 

Taking into account the considerations above, we present an orbit-based procedure for constructing the 
deformation pair of functions. The steps to be followed for deforming two-field interacting models are the 
following: 

1. Choose a deformation function for one of the fields, for example a function fi(fa). Then, as stated in 
([15} we define the deformed field by xi = fi(fa) ( or > equivalently, fa — fx (xi))- For the other field, 
X2 = h{fa) we can write 

fa = [ d 2 fi\ X2 )d X2 = I dx _\ (24) 
J J #2/2 (/ 2 1X2)) 

2. Choose the topological sector to be deformed and an orbit 0(xi,X2) — 0, associated to this sector. Use 
this equation to write Xi as a function of X2, z-e. Xi = F(X2)- 

3. Then impose the condition (|2"2"|) on and use the expression xi = F(X2) to obtain 

dX2 



= ± 



dih [/r 1 (F( X2 ))] 



= ±J [difi\xi)] xl=F[x2) d X 2. (25) 



After integration we obtain fa as a function of X2, which is nothing but the inverse of the deformation 
function f 2 . This is the key result of the present work. 

Thus, the chosen function fi(fa) and the constructed one f 2 (fa) form a pair that takes the original model 
and solutions to a deformed model with a smooth potential, and solutions satisfying the first order equations 
as well as the equations of motion. 

Note that the procedure described above is orbit-dependent and the possibility of finding the second defor- 
mation function is restricted to the ability of explicitly integrating eq. (|25]l . 

V. A DETAILED EXAMPLE 

As an example of application of the extended deformation method we consider the model [2(| [2l| 

W = Xi ~ 3X1 - ^1X2 . re R (26) 
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It presents 4 minima: vaa — [±1,0] on the xi ax is, an d v bb — [0, ±^] ° n the \i ax i s - 
The corresponding first order system of equations is 



dxi _ i _ 2 _ 2 
dx — 1 Xl r A2 



dX2 
da: 



(27) 
-2rxiX2- 



The integrating factor for this system can be found explicitly, and this allows to find the flow-line family of 
curves 

rxl - (1 - 2r)(l - xl) + C X { = 0, r ? ~ (28) 

where C is an integration constant. Real values of C give orbits starting or ending in vaa minima, but there 
exist some critical values. When C = C s = — 2r(±y/r)~ the orbits start and end at different axis, joining 
vaa an d Vbb minima. The existence of these critical values determines several ranges of C in R for which the 
corresponding solutions of (j2"5|) are not kink orbits - see [21| . 

So we have a two-field model and its general orbit equation depending on two parameters (r and C). In 
order to apply the orbit-based procedure described in the previous section, we will consider separately the 
different kind of orbits, corresponding to different regions in parameter space. 



A. Elliptic orbit deformation 



Let us first consider the simplest case in which the integration constant C is taken to be zero and r restricted 
to the interval (0, i). In this case, the orbits are ellipses and (|2"8"]l can be rewritten as 



X2 = F(xi) = ± 



1 - 2r 



(1-Xi) (C = 0, 0<r<i) 



A two-field static solution for the system ([27]) . which satisfies this constraint is [2C 
Xl( x ) — tanh(2rx), 



X2(z) 



1 - 2r 



sech(2rcc) 



(29) 



(30) 



In figure [T] we show the vacua structure and some orbits of the model for different values of the r parameter. 
While the two minima in the horizontal axis form a topological sector (AA-sectov), the two minima in the 
vertical axis (BB-sectoi) cannot be connected by solutions of the first order system (|27l) . 




Figure 1: Minima and orbits of the undeformed model for C = and different r values (r = |, |, j). 
Now, following the prescription established in the preceding section, we choose a deformation function for 



Xi 



0i = A 1 (xi) = arcsinh(xi) 



(31) 



For constructing the deformation function for the other field, \2, we calculate the integral (j2"5|) using 
and CCD 



!>2 = 



dX2 



f{[fiHF( X2 ))] 



2- 



(32) 



27 X 2 



We obtain 



Thus 



02 = / 2 " 1 (X2) 



/ 2 (^2) 



1 - 2r 



■ arcsm 



2(1 - 2r) 



X2 



2(1 - 2r) 
— sin 



1 - 2r 



(33) 



(34) 



With the deformation functions at hand, and making use of (|23[) . we are able to write down the deformed 
potential, which reads 



1 - sinh 2 ^) - 2(1 - 2r) sin 2 (yj^4> 
+ 2 r (1 - 2r) sinh 2 (0! ) tan 2 ( 2 



sech 2 (^i) 



(35) 



Evaluating the fields xi an( i Xi at the solution of the original model we obtain a solitonic solution of the 
deformed model specified by U(<j>i, fo) 



t>i(x) — arcsinh [tanh(2rx)] 



<fo{x) 



1 - 2r . I V2 , .„ , 
arcsm I -— sech(2ra;) 



(36) 



Condition (|22|) is automatically satisfied by this solution. This makes it consistent with the second order 
equations of the deformed system, as can be explicitly verified. 

To show how this deformation acts, in figure [2] we plot both, the deformed and the original solutions. 

X2 (X) ,02 (X) 





-10 -5 5 10 

Figure 2: Solutions of the undeformed (dotted line) and deformed (solid line) models for C = 0, r = |. 
Substituting the original fields by its corresponding deformed partners in (|28[) we obtain the deformed orbit 

cosh 2 (0i ) — 2 cos 2 
which allows writing (f>2 as a function of (j>\ 

= ± \ arccos | cosh(0i ) J + fc7r ( fc € Z ) , (38) 



1 - 2r 



0. 



(37) 
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-6 



Figure 3: Some vacua and orbits of the deformed model with C = 0, r = 1/4. 



and this explicitly shows that the new model presents a periodic vacua structure. Such structure and some 
orbits of the deformed model are shown in fig. [3] 

We can summarize the results obtained for the AA-secioi as follows 



' orb AA 
f 

OrbAA def 
VAA def 

${x) 



r X 2_(i_2r)(l- X 2) =0 ; 
I tanh2rir, \ 1 ~ 2r sech2rx 



v AA = {[l,0],[-l,0]} 



(sinh^y^sin^^)) 



cosh 



2 cos- 



arcsinh (tanh 2rx) , 



= 0. 



(39) 



l-2r 



arcsm 



f ^ sech2ra; 



B. Linear orbits deformation 



It is also possible to find explicit solutions for the BNRT model (|2l)]) in other regions in parameter space. 
For example, for specific values of the integration constant C and parameter r, there are orbits connecting 
one minimun on the xi — ax is with one on the \2 = axis (AS-sectors). As an illustration we address now 
the r = 1 case, and deform orbits for integration constants C — C s — ±2. There are four linear orbits, one 
for each sector, that we label a, 6, c and d, starting from the first quadrant and moving forward clockwise, as 
shown in figure [4] 



X2 




Figure 4: ,45-orbits (C = ±2, r = 1). 



In order to deform these linear sectors we choose again the deformation function for one of the fields as 
Xi = /i(0i) = sinh(</>i), and construct the 72(^2) by using the corresponding orbits. The resulting vacua 
structure and orbits for the deformed models are shown below. 



For the a sector we obtain 



f orb^l 
X {a \x) 
/(-) 

° rh( ABdef 

(a) 
V AB def 



X2+Xi-1 = 0; ^ = {[1,0], [0,1]} 

( |(1 - tanha;), |(1 + tanha;) ) 
( sinh + sinh (f>2 ) 

sinh <p2 + sinh <f>\ = 

{[0,0],[0,-sinh- 1 (2)],[sinh- 1 (l),sinh- 1 (-l)],[sinh- 1 (-l),sinh- 1 (-l)]} 
( — sinh -1 1(— 1 + tanha;), sinh" 1 | (—1 + tanha:) ) 




Xi _! 



(a)Original orbit (b)Deformed orbit 

Figure 5: Deformation of AB {a) sector (C = -2, r = 1). 




1 / 



Figure 6: Compared profiles of original (dashed line) and deformed (solid line) defect solutions. AB^ 



For the b sector we get 



( orb% 

fib) 

0rb( ABdef 



(b) 
' ] AB def 

I $ {b) {x) 



X2-Xi + 1 = 0; v%= {[1,0], [0,-1]} 

( ^(1 — tanha;), — 1(1 + tanha;) ) 
( sinh 01 j — 1 + sinh 2 ) 

sinh 4>2 — sinh <f)\ = 

{[0,0], [O^inlT 1 2], [sinh _1 (l),smh _1 (l)], [sinh _1 (-l), sinh _1 (l)]} 
( — sinh 1 i(— 1 + tanha;), — sinh" 1 s(— 1 + tanha;) ) 



X2 
It 




Xl 



1 

• 




-1 


-1 




l 



-1 

(a)Original orbit (b)Deformed orbit 

Figure 7: Deformation of AB (b) sector (C = 2, r = 1). 



01 



xj b) (x), 0i b| (x) 




Figure 8: Compared profiles of original (dashed) and deformed (solid) defect solutions. AB^ sector 



the case of the c sector we obtain 



( orb% 

/(c) 



orfe 



ABdef 



AB def 
I (C) (^) 



«il = {[-1,0], [0,-1]} 



X2 + Xi + 1 = 0; u ab 
( — 1(1 — tanhz), — 1(1 + tanhx), ) 
( sinh </>! , — 1 + sinh </> 2 ) 

sinh (f>2 + sinh 0i = 

{[0,0], [0,sinh _1 2], [sinh _1 (l), sinh _1 (l)], [sinrT^-l), sinh _1 (l)]} 
( sinh -1 | (— 1 + tanhcc) , — sinh -1 ~ (— 1 + tanhrr) ) 



X2 

it 



02 




Xl 




(a)Original orbit. 



(b)Deformed orbit. 



Figure 9: Deformation of AB (c) sector (C = 2, r = 1). 
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X{ c) (x) , 0< c) (x) 



-5 



-1 




Figure 10: Compared profiles of original (dashed line) and deformed (solid line) defect solutions. AB^ sector. 



Finally, for the d sector we have 



( orb% 

(d) 



■■ X2-Xi-1 = 0; vf B = {{-!,% [0,1]} 

X {d \x) = (-±(1 - taimar), |(1 + tanha;),) 



0rh{ ABdef 1 

v [d) 

AB def — 



( sinh 0i,l + sinh <p2 
sinh <f>2 — sinh = 

{[0,0], [0,-sinlT 1 2], [sinh _1 (l),sinh -1 (-l))], [sinh _1 (-l), sinh _1 (-l)]} 
(sinh^ 1 | (—1 + tanhs) , sinh^ 1 | (—1 + tanha;)) 



(43) 




(a)Original orbit. (b)Dcformed orbit. 

Figure 11: Deformation of AB {i) sector (C = —2, r = 1). 
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Xl d) (x), 0< d) (x) 



2"/ 



X j d) (x), 0« d »(x) 



t x 




1 

1 

~2~/ 






-5 










1 








~ 2/ 





Figure 12: Compared profiles of original (dashed line) and deformed (solid line) defects. AB^> sector. 



In the above examples, we have considered two cases in which C has been fixed. Let us now move to other 
regions in parameter space, letting C to be undetermined. This will of course increment the complexity of the 
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problem but, as we will show below, there is still a lot of possibilities of generating new models. This illustrates 
the richness of results that can be obtained by exploiting the deformation method extended to interacting field 
models. 



C. An integrable case 



Recall the general orbit equation 



x? = i- 



C i 
■Xi- 



l-2r z l-2r 

Now, following the prescription for generating the deformed model, we are led to 

Xi = fiifa) = sinh 0i ; fa = /f^Xi) = arcsinhxi 

and 

dX2 



2 l-2rX2 l-2r ^2 



(44) 



(45) 



(46) 



Integral (|46|) is highly non trivial for arbitrary values of C and r. However it is integrable for some values 
of r, in particular for r = \- For this value there is only one critical constant C s = — and kink orbits arise 
if C € [C s ,oo). In this case, ([46]) is an elliptic integral, that can be written in the form 



dX2 



\ /Z 2C J ^( X 2 2 -Xi)(xl~^y 



(47) 



with Ai = —± (1 + v 7 ! + 64C) and A 2 = (l - y/l + 64C) . 

By an appropriate change of the integration variable, integral (|47[) can be solved in terms of Jacobian elliptic 
functions [22|. The result is 



1 



Aiv^C 



X.2 



(48) 



We note that for all the Jacobi elliptic functions appearing in this work we have to take its real part, as we 
are dealing with real scalar fields and they are solutions of a physical problem. 

This solution presents distinguishable behaviors depending on C taking values on the regions (— g^,— 
(— 0) or (0, oo). For the special values C = and C = — the deformed system presents field solutions in 
terms of elementary functions rather than Jacobi elliptic functions. 

Putting the original field solutions in terms of the deformed ones 



Xi = sinh fa , X2 = V Aa sn V AiV-2C fa 



(49) 



we can write the explicit form the deformed model potential. We obtain 



U(fa,fa;\, C) = \ 



1 - sinli 



fa - A 2 sn 2 (v/aTv^C fa- ^A 2 /Ai) 



sech 



■ — - — sinh 2 fa 



scr 



/A7V 3 2C0 2 ; v / VAi') 



dn 2 



(50) 



A field solution of the undeformed system ([27]) , satisfying the orbit |44|) for r = 4 is 



(51) 
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where b 2 = l/yl + 32C. Therefore, using the deformation functions (|4"5|) . the solutions for the deformed 
system read 



bi (x) = arcsinh 



sinh(x) 



b 2 + cosh(x) 



1 



V-2AiC 



b/V^2 



y/cosh{x) + b 2 ' V -^1 



(52) 



In figure [TBI we plot the original and deformed profiles of the field solution (taking C = 1). It is remarkable 
how similar is the behavior of these solutions and the ones plotted in figure despite the much more involved 
analytical expressions in this last case (compare formula (|50[) with formula ([35])). 



xf =1 (x), ^ =1 (x) 





Figure 13: Deformed (solid) and undeformed (dotted) solutions for BNRT model (r = |, C = 1). 



From the results obtained above we stress that the orbit-based deformations can be applied to every kink 
orbit of the BNRT model (|26|) for the integrable case of r — 1/4. For almost any value of the constant C 
we obtain different models in terms of elliptic Jacobi functions, with two exceptions: C — and C — —1/64. 
Also, we note that there is only one kink solution for each member of the family of deformed models, coming 
from the appropriate kink orbit in the original model. 



D. Another integrable case. 



For r — 1 the deformed model for every kink orbit can also be found, given that, as for the former case, it 
is possible to integrate (|46|) for any value of C. In fact, setting r = 1 in (|46| and taking <p\ — arcsinh(xi), we 
obtain the deformed cf>2 field 



1 



xi 



(53) 



Writing the original field solutions in terms of the deformed ones, we can find the explicit form for the 
deformed model potential, which turns out to be 

2 ) 2 ] 2 

(54) 



[/(0 1 ,0 2 ;r = l,C) 



\ sech(0i) 2 [l - sinh(^i) 

8+4. - 



+isinh(0 1 ) 2 e- 2 ^ 
The solutions of the BNRT model for r 

XX 



64 

* -4 C e*^ 2 +C 2 



e "'- 2 (-8 + 4e 2<fc -4Ce* a 



8e-*2+4e*2-C 2 

1 can be written as 



(e 2 *) 2 - C 2 + C 2 



(e 2 * - Ci) 2 - CV 



X2 



2C 2 e 2x 



{e 2x -df 



CI 



(55) 



with C2 = 2C\/C and |C| > \C S \ = 2. As was said before, the solutions given in formula (|55p are not kinks 
if \C\ < 2, as both \i and X2 become singular for this values of C. 
The corresponding deformed solutions read 



arcsinh 



((i^-cl+cl\ 



(56) 
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For the specific value C = 2, the solutions ([55]) are both kink-like, as well as the deformed solutions (|56| . 
and its profiles are very similar to that of the case C = — 2, plotted in Fig. [TUl 

For C > 2 we find that the profiles of xi and X2 in (|55p are kink and lump-like respectively, as shown in 
Fig. M 



Xf- 2 (x) , <t>f 2 (x) X f 2 (x) , 0f 2 (x) 




Figure 14: Deformed (solid) and original (dotted) solutions for BNRT model for r = 1, Ci = —1 and C = 3. 



Before ending this section, let us recall that the values r = 1/4 and r = 1, for which the integral (|46|) can 
be analyticaly implemented, present special features: kinks of the r = 1/4 family are separatrix trajectories 
of a mechanical problem which is Hamilton- Jacobi separable in parabolic coordinates, and for r = 1, the 
mechanical problem is separable in Cartesian coordinates; these problems are known as Liouville type III and 
IV, respectivelly (23|. However, even though the separable Liouville type problems seem to lead to integrable 
expressions in equation (|25p . in practice not all the kink orbits are algebraic, and this may preclude the presence 
of analytical solutions. 



VI. COMMENTS AND CONCLUSIONS 

In this work we have presented a generalization of the deformation method, first introduced in Q, which 
allows to generate deformed potentials and solutions given a model of two real scalar field and its solutions. 
The main new ingredient consist in the need of imposing a constraint on the functions used to deform the 
fields, required to preserve the relation between the original solutions, that live in orbits of the configuration 
space. As the construction of the deformed solutions involves the orbit constraint, the deformation must be 
implemented independently on each topological sector. Consequently, different orbits (even when belonging 
to the same sector) can lead to different full deformed models and solutions. 

The present version of the method applies to models with an associated superpotential W{\i, Xi)- Although 
we have considered deformation functions depending on a single field, more general functions can be used, and 
reduced to the former case with the use of the orbit itself. Deformation functions depending on both scalar 
fields are now being considered in the context of models with holomorphic superpotentials f24j , as is the case 
of the bosonic sector of the (l+l)-dimensional N=2 SUSY Wess-Zumino model. We will further report on this 
possibility in future work. 

As the described procedure is orbit dependent, whenever deforming a two-field model, the integrability of 
equation (f2"5|) has to be analyzed separately in each case. However, in case of two-field theories which are 
associated to separable mechanical systems of the Liouville type, it seems that the present method will work 
very nicely, at least when one is restricted to consider algebraic kink orbits. Interesting examples of this kind 
are known, as the celebrated MSTB model [25| and other models proposed recently in Ref. [261 ]. 

This new version of the deformation procedure provides a tool for studying more sophisticated systems. It 
allows one to generate a diversity of new systems with their corresponding solutions, which may contribute to 
improve the understanding of complex problems. 
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